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Abstract 



In this article we present an analytical method for deriving the relationship be- 
tween the pressure drop and flow rate in laminar flow regimes, and apply it to the 
flow of power-law fluids through axially-symmetric corrugated tubes. The method, 
which is general with regards to fluid and tube shape within certain restrictions, 
can also be used as a foundation for numerical integration where analytical ex- 
pressions are hard to obtain due to mathematical or practical complexities. Five 
converging-diverging geometries are used as examples to illustrate the application 
of this method. 



1 INTRODUCTION 
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1 Introduction 

Modeling the flow in corrugated tubes is required for a number of technological and 
industrial applications. Moreover, it is a necessary condition for correct description 
of several phenomena such as yield-stress, viscoelasticity and the flow of Newtonian 
and non- Newtonian fluids through porous media [1-4]. In the literature of fluid 
mechanics there are numerous studies on the flow through tubes or channels of 
corrugated nature. Most of these studies use numerical techniques such as spectral 
and finite difference methods [1, 5-11]. Some others adopt analytical approaches 
based on simplified assumptions and normally deal with very special cases [12, 13]. 
The current paper presents a mathematical method for deriving analytical relations 
between the volumetric flow rate and pressure drop in corrugated circular capillaries 
of variable cross section, such as those illustrated schematically in Figure 1. It 
also presents five examples in which this method is used to derive flow equations 
for non-Newtonian fluids of power-law rheology. The method, however, is more 
general and can be used for other tube shapes and other rheologies. The following 
derivations assume a laminar flow of a purely- viscous incompressible fluid where the 
tube corrugation is smooth and relatively small to avoid complex flow phenomena 
which are not accounted for in the underlying assumptions of this method. 




Figure 1: Profiles of converging-diverging axisymmetric capillaries. 



2 P-Q RELATION FOR POWER-LAW FLUIDS 
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2 P — Q Relation for Power-Law Fluids 

The power-law, or Ostwald-de Waele, is a widely-used fluid model to describe the 
rheology of shear thinning fluids. It is one of the simplest time-independent non- 
Newtonian models as it contains two parameters only. The model is given by the 
relation [14, 15] 

= - = cr-' (1) 

7 

where /z is the fluid viscosity, r is the stress, 7 is the strain rate, C is the consistency 
factor, and n is the flow behavior index. In Figure 2 the bulk rheology of this model 
for shear-thinning case is presented generically on log-log scales. Although the 
power-law is used to model shear-thinning fluids, it can also be used for modeling 
shear-thickening by making n greater than unity. The major weakness of power- 
law model is the absence of plateaux at low and high strain rates. Consequently, 
it fails to produce sensible results in these flow regimes. 

A 




' ► 

Figure 2: The bulk rheology of power-law fluids on logarithmic scales for flnite 
strain rates (7 > 0). 



2.1 Conic Tube 
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The volumetric flow rate, Q, of a power law fluid through a circular capillary of 
constant radius r with length x across which a pressure drop P is applied is given 
by 

^ ~ 2xCi/«(3n + 1) \V) 
A derivation of this relation, which can be obtained from Equation 1, can be found 
in [16, 17]. On solving this equation for P the following relation is obtained 

^ ^ 2CQ^{3n + irx 

For an infinitesimal length, 5x, of a capillary, the infinitesimal pressure drop for 
a given flow rate Q is 

^ 2Cg"(3n + irfe 

For an incompressible fluid, the volumetric flow rate across an arbitrary cross 
section of the capillary is constant. Therefore, the total pressure drop across a 
capillary of length L with circular cross section of varying radius, r(x), is given by 

p ^ 2Cg"(3n + l)" /-^ 

In the following sections. Equation 5 will be used to derive analytical expressions 
for the relation between pressure drop and volumetric flow rate for flve converging- 
diverging geometries of axially-symmetric pipes. The method can be equally ap- 
plied to other geometries as long as the tube radius can be expressed analytically 
as a function of the tube axial distance. 



2.1 Conic Tube 
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Figure 3: Schematic representation of the radius of a conically shaped converging- 
diverging capillary as a function of the distance along the tube axis. 

2.1 Conic Tube 

For a corrugated tube of conic shape, depicted in Figure 3, the radius r as a function 
of the axial coordinate x in the designated frame is given by 



r{x) = a + b\x\ -L/2<x<L/2 a,b>0 



(6) 



where 



a = Rr. 



and 



(7) 



Hence, Equation 5 becomes 



2CQ"(3n+ 1)" 



L/2 



dx 



-L/2 



2CQ"(3n + 1)" 

4Cg"(3n + 1)" 
37r"n"+i6 



(a + 



3bn{a — bx)^" 
1 



+ 



-L/2 



(a + 6L/2)3"_ 



3bn{a + 



(9) 
(10) 



that is 



P 



2LCQ"(3n + 1)" 

37l^n"'^^(^Rmax Rn 



pSn D3n 
min ''max 



2.2 Parabolic Tube 
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2.2 Parabolic Tube 
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Figure 4: Schematic representation of the radius of a converging- diverging capillary 
with a parabolic profile as a function of the distance along the tube axis. 



For a tube of parabolic profile, depicted in Figure 4, the radius is given by 



r{x) = a + hx^ 



L/2<x<L/2 a,b>0 



(12) 



where 



and 



L 



(13) 



Therefore, Equation 5 becomes 



2CQ"(3n + r 



L/2 



dx 



2CQ"(3n + 1) 



-L/2 (a + hx'^f'"'^^ 

3n 



^'^Tfe;^ ■^^■2-3" + !' 2'-- 



L/2 



-L/2 



(14) 



(15) 



where 2-^1 is the hypergeometric function. Therefore 



^ 2LCQ" 3n + l " ^ /I 3 

P = . ^' 2i^i -, 3n + 1; -; 1 — 



(16) 



2.3 Hyperbolic Tube 
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2.3 Hyperbolic Tube 

For a tube of hyperbolic profile, similar to the profile in Figure 4, the radius is 
given by 



r{x) = Va + bx^ - L/2 <x < L/2 a,b>0 



[17) 



where 



p2 



and 



L 



Therefore, Equation 5 becomes 



P 



2CQ"(3n + 1)" 



2CQ"(3n + 1) 



L/2 



dx 



L/2 (a + hx 
^ + 1 



2N(3n+l)/2 



3n + l 
2 



(19) 



X 



a + hx"^ 



1 3n + 1 3 hx' 



T 2 '2' a 



L/2 



(20) 



-L/2 



where 2-^1 is the hypergeometric function. Therefore 



^ ^ 2LCQ"(3n + l)" (^\^ 3n + 1 _ 3 _ ^ 



7r"n«i?i"+^ ' ^2' 2 '2' 



. 

mm 



(21) 



2.4 Hyperbolic Cosine Tube 

For a tube of hyperbolic cosine profile, similar to the profile in Figure 4, the radius 
is given by 



r(a;) = acosh(6s) - L/2 < x < L/2 a>0 



(22) 



where 



and 



2 1 -Rrr 

— arccosh I — 

L V iXr, 



(23) 



2.5 Sinusoidal Tube 
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Hence, Equation 5 becomes 



P 



2Cg"(3n + l)'^ /"^/^ dx 



7r"n"a3"+i J^^/^ cosh^"+^(6, 



X) 



2CQ''{3n + 1)" 

3^n^n+lQ3n+15 



sinh(6x) 



cosh'^" (bx)\/ — sinh^ (6a; 



1 3n 2 — 3n , 2/7 ^ 
1 2' — 2"' — 2 — ' ' 



On evaluating the last expression and taking its real part which is physically 
significant, the following relation is obtained 



2LCQ"(3n + 1) 



37rn^n+i^^^^^3n^arccosh (f^ 



\ 2 L \ , , , ^2 



(26) 



where Im(2-Fi) is the imaginary part of the hypergeometric function. 



2.5 Sinusoidal Tube 




Figure 5: Schematic representation of the radius of a converging- diverging capillary 
with a sinusoidal profile as a function of the distance along the tube axis. 



For a tube of sinusoidal profile, depicted in Figure 5, where the tube length L 
spans one complete wavelength, the radius is given by 



r{x) = a-bcos{kx) - L/2 < x < L/2 a>b>0 (27) 



2.5 Sinusoidal Tube 
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where 



^^Rrna^^rran ^ ^28) 



Hence, Equation 5 becomes 



2Cg«(3n + 1)' 



n fL/2 



dx 



vr"n" y_^/2 [a-h cos {kx)] 



2CQ"(3n + 1)" 



3n+l 



(29) 



sin^(A;x) 



sin(/cx) [a — 6 cos (A;x)] " y — IP' 



21 



L/2 



(30) 



-L/2 



where 21 is the Appell hypergeometric function given by 



1 1 



a — h cos{kx) a — b cos{kx) 



2l = Fi -3n;-,-;l-3n; 

22 a + o a — 



(31) 



On taking the hmit and X — )■ ^ and considering the real part, 

the following relation is obtained 



2L(:7Q"(3n + 1)" 



u/i 11 -I V -"'maa::-' i-ri 



1 1 



=Im( Fi ( -3n;-,-;l-3n;l, 



(32) 



where Im(Fi) is the imaginary part of the Appell hypergeometric function. 



It is noteworthy that all these relations (i.e. Equations 11, 16, 21, 26 and 32), 
are dimensionally consistent. Moreover, they have been thoroughly tested and 
validated by numerical integration. 



3 CONCLUSIONS 
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3 Conclusions 

In this article a mathematical method for obtaining analytical relations between 
the flow rate and pressure drop in axisymmetric corrugated tubes is presented and 
applied to the flow of power-law fluids. The method is illustrated by flve examples 
of circular capillaries with converging-diverging shape. This method can be used 
to derive similar relations for other types of fluid and other types of tube geometry. 
It can also be used as a base for numerical integration where analytical relations 
are difficult to obtain due to mathematical or practical reasons. 



3 CONCLUSIONS 



Nomenclature 

7 strain rate (s~^) 

II fluid viscosity (Pa.s) 

T stress (Pa) 

C consistency factor (Pa.s") 

Fi Appell hypergeometric function 

2F1 hypergeometric function 

L tube length (m) 

n flow behavior index 

P pressure drop (Pa) 

Q volumetric flow rate (m^.s~^) 

r tube radius (m) 

Rmax maximum radius of corrugated tube (m) 

Rmin minimum radius of corrugated tube (m) 

X axial coordinate (m) 
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